Dysregulation of microRNAs and tRNA-derived ncRNAs in mesothelial and mesothelioma cell lines after asbestiform fiber exposure

Experimental evidence demonstrated that fluoro-edenite (FE) can develop chronic respiratory diseases and elicit carcinogenic effects. Environmental exposure to FE fibers is correlated with malignant pleural mesothelioma (MPM). An early diagnosis of MPM, and a comprehensive health monitoring of the patients exposed to FE fibers are two clinical issues that may be solved by the identification of specific biomarkers. We reported the microRNA (miRNA) and transfer RNA-derived non coding RNA (tRNA-derived ncRNA) transcriptome in human normal mesothelial and malignant mesothelioma cell lines exposed or not exposed to several concentration FE fibers. Furthermore, an interactive mesothelioma-based network was derived by using NetME tool. In untreated condition, the expression of miRNAs and tRNA-derived ncRNAs in tumor cells was significantly different with respect to non-tumor samples. Moreover, interesting and significant changes were found after the exposure of both cells lines to FE fibers. The network-based pathway analysis showed several signaling and metabolic pathways potentially involved in the pathogenesis of MPM. From papers analyzed by NetME, it is clear that many miRNAs can positively or negatively influence various pathways involved in MPM. For the first time, the analysis of tRNA-derived ncRNAs molecules in the context of mesothelioma has been made by using in vitro systems. Further studies will be designed to test and validate their diagnostic potential in high-risk individuals' liquid biopsies.

www.nature.com/scientificreports/ Analysis (PCA) showed a separation between MeT-5A and JU77 cell lines in their miRNAs and tRNA-derived ncRNAs expression. The greater separation was observed between the two different cell lines (PC1: 45.43%) while treated and untreated cell lines showed minor differences (PC2: 18.16%) ( Fig. 2A). Moreover, to better highlight intra-cell line differences in small ncRNA expression, we performed PCAs specifically for each cell line (Figs. 2B, C). These confirmed a great separation between treated and untreated JU77 (PC1: 44.95%) (Fig. 2B).
As about MeT-5A cells, the separation was greater by comparing untreated cells with that one exposed to 50 µg/ ul of FE fibers (PC1: 29.22%) (Fig. 2C). In combination, these results seem to suggest that although differences in small ncRNA expression between treated and untreated MeT-5A were detected, major effects were observed in the JU77 cell line. The differentially expressed miRNAs and tRNA-derived ncRNAs were: 1171 in untreated JU77 vs. MeT-5A, 960 in JU77 vs. MeT-5A exposed to 10 ug/ml FE fibers, and 969 in JU77 vs. MeT-5A exposed to 50 ug/ml FE fibers (Fig. 3A). The common population of differentially expressed miRNAs and tRNA-derived ncRNAs between the two cell lines increased with the exposure to FE fibers.
When we compared miRNAs and tRNA-derived ncRNAs expression between unexposed vs. exposed MeT-5A the results showed several differentially expressed molecules. In particular, miRNAs and tRNA-derived ncRNAs differentially expressed in MeT-5A untreated vs. exposed to 10 and 50 ug/ml FE fibers were 7 and 64, respectively (Fig. 3B). Indeed, the miRNAs and tRNA-derived ncRNAs differentially expressed in MeT-5A increased in dose-dependent manner with the exposure to FE fibers. On the other hand, the comparison between untreated vs. FE-treated JU77 showed a substantial difference between the differentially expressed miRNAs and tRNAderived ncRNAs. In particular, miRNAs and tRNA-derived ncRNAs differentially expressed in JU77 unexposed vs. exposed to 10 and 50 ug/ml FE fibers were 549 and 579, respectively (Fig. 3B).
In more detail, the down-regulated miRNAs and tRNA-derived ncRNAs were: 6 and 33 in untreated MeT-5A vs. exposed to 10 and 50 ug/ml FE fibers, respectively (Fig. 4C). While in the case of neoplastic cells, the downregulated miRNAs and tRNA-derived ncRNAs were: 183 and 186 in untreated JU77 vs. exposed to 10 and 50 ug/ JU77 controls vs. exposed to 10 and 50 ug/ml FE fibers; (C) MeT-5A controls vs. exposed to 10 and 50 ug/ml FE fibers.   4C). Among these, 3 were in common between the MeT-5A untreated and exposed to 10 ug/ml FE fibers vs. MeT-5A untreated and exposed to 50 ug/ml FE fibers (tRF-42-YDU37Y93Y086Z34J2, tRF-40-XENDBP1IUUK7VZ0R, and tRF-41-XENDBP1IUUK7VZ0RB), 113 were in common between the JU77 untreated and exposed to 10 ug/ml FE fibers vs. JU77 untreated and exposed to 50 ug/ml FE fibers (87 tRNAderived ncRNAs and 26 miRNAs), only hsa-miR-3618 was in common between the MeT-5A untreated and exposed to 50 ug/ml FE fibers vs. JU77 untreated and exposed to 10 ug/ml FE fibers, no compound was in common between the MeT-5A untreated and exposed to 50 ug/ml FE fibers vs. JU77 untreated and exposed to 50 ug/ ml FE fibers (Fig. 4C, Table 4). If we compare MeT-5A exposed to 10 ug/ml FE fibers vs. MeT5A untreated and MeT-5A exposed to 50 ug/ml FE fibers vs. MeT5A untreated and JU77 exposed to 10 ug/ml FE fibers vs. JU77 untreated the only compound in common was hsa-miR-1248 (Fig. 4C, Table 4). tRF-40-8B7OIQ1IBM4LV7U6, tRF-31-79MP9P9MH57SD, tRF-32-P4R8YP9LON4V3, tRF-25-SP58309MUK were the down-regulated tRNAderived ncRNAs in common among MeT-5A exposed to 50 ug/ml FE fibers vs. MeT5A untreated and JU77 exposed to 10 ug/ml FE fibers vs. JU77 untreated and JU77 exposed to 50 ug/ml FE fibers vs. JU77 untreated (Fig. 4C, Table 4). Among all samples, tRF-33-79MP9P9MH57SD3 was the only down-regulated compound in common (Fig. 4C, Table 4). As for the up-regulated miRNAs and tRNA-derived ncRNAs, these were: 1 and 31 in untreated MeT-5A vs. exposed to 10 and 50 ug/ml FE fibers, respectively (Fig. 4D). While in the case of neoplastic cells, the up-regulated miRNAs and tRNA-derived ncRNAs were: 366 and 393 in untreated JU77 vs. exposed to 10 and 50 ug/ml FE fibers, respectively (Fig. 4D). Among these, only ts-96 was in common between the MeT-5A untreated and exposed to 10 ug/ml FE fibers vs. MeT-5A untreated and exposed to 50 ug/ml FE fibers, 289 were in common between the JU77 untreated and exposed to 10 ug/ml FE fibers vs. JU77 untreated and exposed to 50 ug/ml FE fibers (among these only one was a miRNA, specifically hsa-miR-615-3p) (Fig. 4D, Table 5). tRF-16-PSSELQB, and tRF-21-JYSWRYVMD. MeT-5A cell lines both not exposed and exposed to FE fibers; (C) The same cellular line both not exposed and exposed to FE fibers. Venn diagrams showing up-regulated miRNAs and tRNA-derived ncRNAs distribution between: (B) JU77 vs. MeT-5A cell lines both not exposed and exposed to FE fibers; (D) The same cellular line both not exposed and exposed to FE fibers. www.nature.com/scientificreports/ Table 2. Down-regulated ncRNAs in common among JU77 (human malignant mesothelioma) vs. MeT-5A (human normal mesothelium) controls and JU77 vs. exposed to 10 and 50 ug/ml FE fibers; down-regulated ncRNAs in common among JU77 vs. MeT-5A exposed to 10 and exposed to 50 ug/ml FE fibers.
Surely, between miRNAs and tRNA-derived ncRNAs molecules, the latter were found to be the most deregulated among the samples compared.
The results of the differential expression analysis were reported in Sup. Table 1.

Pathways analysis.
Once we identified the differentially expressed miRNAs for each comparison, we tried to investigate the impact of their dysregulation in metabolic and signaling pathways by using MITHrIL 50 .
MITHrIL fully exploits the topological information encoded by pathways when computing perturbation scores. Pathways are then modeled as complex graphs where each node is a biological element (protein-coding gene, miRNA, or metabolite), and each edge is an interaction between them 50 . Importantly, MITHrIL takes into account experimentally validated miRNA-mRNA interaction in order to predict their effects on biological pathways 50 .
The results demonstrated clear patterns of negative and positive perturbation scores involving 39 different pathways. Considering untreated JU77 vs. MeT-5A, 24 pathways showed a negative perturbation score while 15 pathways showed a positive perturbation score.
A similar trend was found in the case of the JU77 vs. MeT-5A exposed to 10 ug/ml FE fibers, indeed 23 pathways showed a negative perturbation score while 16 pathways showed a positive perturbation score. In fact, the Jak-STAT signaling pathway, the toll-like receptor signaling pathway, and the thyroid hormone signaling pathway showed a different trend between the samples mentioned above. In particular, in samples not exposed and exposed to 10 ug/ml FE fibers, the Jak-STAT signaling pathway and the thyroid hormone signaling pathway showed a negative and positive perturbation score, respectively. Both of these pathways reverse their score again in the samples exposed to the FE fiber concentration equal to 50 ug/ml. On the contrary, the toll-like receptor signaling pathway showed a positive perturbation score in the samples not exposed and a negative perturbation score in the samples exposed to both FE fibers concentrations.
Considering the JU77 vs. MeT5A exposed to 50 ug/ml FE fibers, 26 pathways showed a negative perturbation score while 13 pathways showed a positive perturbation score. Comparing these results with the cells that did not undergo any exposure to FE fibers, the phospholipase D signaling pathway and the toll-like receptor signaling pathway showed a negative perturbation score after FE fibers exposure (Fig. 5).
The pathway analysis performed between untreated vs. FE fibers treated JU77 and MeT-5A and the main impacted pathways showed clear patterns of positive correlations involving 23 different pathways. Only the NFkappa B signaling pathway showed negative perturbation in the neoplastic samples (Fig. 6).
The results of the pathways analysis were reported in Sup. Table 2. Table 3. Up-regulated ncRNAs in common among JU77 (human malignant mesothelioma) vs. MeT-5A (human normal mesothelium) controls and JU77 vs. exposed to 10 and 50 ug/ml FE fibers; down-regulated ncRNAs in common among JU77 vs. MeT-5A exposed to 10 and exposed to 50 ug/ml FE fibers.  Fig. 7. Therefore, with the exception of the presence of fibers, neoplastic cells are not distinguished from healthy cells by cellular morphology. In accordance with Panzetta et al. 51 morphology alone is not sufficient to discriminate malignant cells from benign cells. In addition, the in vitro viability of MeT-5A and JU77 from 6 to 72 h exposure to FE has been evaluated to examine the effects of FE fibers upon mesothelium and MM. The results showed that MeT-5A cells were more sensitive to FE fibers compared to JU77 cells, as shown in Fig. 8.
Mesothelioma knowledge-network inference. NetME 47 has been used to recover known interactions among malignant mesothelioma, asbestos, miRNAs, and other biological elements. For this purpose, we selected     Table 2 shows the most relevant sentences from 38 papers analyzed by NetME to infer the knowledge network. From these documents, it is clear how several miRNAs were involved in apoptosis, necrosis, and cell growth inhibition. Among these, the main miRNAs studied are miR-302b, miR-192, and miR-193a. The miR-193b is involved in the autophagy process while miR-1, miR-34a, miR-215, miR-16, miR-320a, and miR-21 seem to participate as oncogenes or oncosuppressors in malignant mesothelioma. To date, in literature, there are no papers that have analyzed tRNA-derived ncRNAs in malignant mesothelioma both induced and not induced by FE fibers. Data collected by NetME showed several research papers that demonstrated the correlation between asbestos exposure and different cancers. In fact, breathing asbestos fibers causes not only malignant mesothelioma, but also may lead to laryngeal cancer, prostate cancer, colorectal cancer, lung cancer, and ovarian cancer. Malignant mesothelioma and lung cancer have also been associated with exposure to carbon nanotubes. In Fig. 9, is shown the mesothelioma-miRNAs interactions network achieved by running the following cypher query against the whole mesothelioma network Each path has been inferred through the NetME syntactic analysis performed over the mesothelioma-based documents collected in PubMed Central.

Discussion
Small RNA-Seq transcriptome profiling of healthy mesothelium and MPM in vitro has been evaluated to highlight the deregulated miRNAs and tRNA-derived ncRNAs and the different pathways involved in an aggressive cancer such as MPM. The results showed that MeT-5A cells were more sensitive to FE fibers compared to JU77 cells and that morphology alone is not sufficient to discriminate malignant cells from benign cells. Certainly, a big difference was found in the number of deregulated miRNAs and tRNA-derived ncRNAs between tumor and non-tumor samples both exposed and not exposed to FE fibers. The common population of differentially expressed miRNAs and tRNA-derived ncRNAs between the two cell lines increased with the exposure to FE fibers. Though, the effect of exposure to FE fibers is more evident in the expression of miRNA and tRNA-derived ncRNAs in tumor samples than in non-tumor samples. Within the same cell line, the effect of exposure to FE fibers caused a dose-dependent increase in differentially expressed miRNAs and tRNA-derived ncRNAs. While the increase was noticeable in MeT-5A cells, in JU77 cells this was less pronounced. In the latter, in fact, the differentially expressed miRNAs and tRNA-derived ncRNAs were already many even in the neoplastic cells that had not undergone treatment with the cancer-causing fibers.
Several molecules aroused interest in their behavior in the various samples analyzed. Among them is a list of down-regulated and up-regulated miRNAs and tRNA-derived ncRNAs in common both between the two cell lines, and the different conditions of exposure with the FE fibers. While tRNA-derived ncRNAs have never been explored in the MPM, miRNAs have already been studied in this context. Several miRNAs are considered potential tumor suppressors in MPM. Expression of miR-15/16 family was consistently down-regulated both in MPM specimens and in cell lines and had tumor suppressor function in MPM 58 , in accordance with our results. miR-193a-3p is also considered a potential tumor suppressor in MPM 59 , indeed these miRNAs belonging to the group of mesomiRs (MM-associated miRNAs) 60,61 could be considered a novel therapeutic approach for MPM.
In the present paper, we analyzed several pathways that are involved in the pathogenesis of malignant mesothelioma. It is very interesting to point out the strong perturbation scores involving the above-mentioned pathways in MPM vs. healthy mesothelial cells. Among these, there was the involvement of pathways that have important functions in inflammatory processes and in angiogenesis. Pathways that showed a different trend in healthy vs tumor samples at different FE fiber exposures were Jak-STAT, phospholipase D, toll like receptors, and www.nature.com/scientificreports/ thyroid hormone signaling pathways. Jak-STAT signaling pathway is involved in processes such as cell death, and carcinogenesis 62 . Phospholipase D activity is significantly increased in cancer tissues and cells, indicating that it plays a critical role in signal transduction, cell proliferation, and anti-apoptotic processes. In addition, phospholipase D is a downstream transcriptional target of proteins that contribute to inflammation and carcinogenesis such as Sp1, NF-kappa B, TCF4, ATF-2, NFATc2, and EWS-Fli 63 . Instead, toll like receptors have roles  www.nature.com/scientificreports/ and functions in anti-cancer immunity and in tumor rejection 64 . Thyroid hormones not only regulate the physiological processes of normal cells but also stimulate cancer cell proliferation via dysregulation of molecular and signaling pathways 65 . Inflammation plays a central role since mesothelioma is a multicentric neoplasm, which originates from inflammatory foci. Inflammation has been correlated with cancer, enhancing the development of malignancies 66 . In particular, the chemokine and TGF-beta signaling pathways lead to acute and chronic inflammation, the latter resulting in several fiber-associated pulmonary and pleural diseases 14,67 . The inflammasome is responsible for the activation of inflammatory processes via multiple mechanisms 67 that induce a process of cell death called pyroptosis, characterized by both apoptosis and necrosis. Cell death mechanisms and release of chemokines and cytokines may help cancer regress, and resist toxicity by fibers and cell growth 68 , in inflammasome-dependent and -independent pathways 69 . Indeed, apoptosis is a mechanism for removing cells with irreparable damage induced by FE fibers, causing genetic changes that predispose cells to a neoplastic transformation 70 . Specific signaling pathways such as sphingolipid, FoxO, and Hippo pathways are involved in many important signal transduction processes such as cell proliferation, and apoptosis [71][72][73] . Dysregulation of the Hippo signaling pathway is highly conserved by phosphorylating and inhibiting the transcription co-activators YAP and TAZ, key regulators of proliferation and apoptosis. On the contrary, dephosphorylated YAP/TAZ translocates into the nucleus and activates gene transcription through binding to TEAD family and other transcription factors. Such changes in gene expression promote cell proliferation and stem cell/progenitor cell self-renewal but inhibit apoptosis, thereby promoting tissue regeneration, and tumorigenesis 72 . An experimental model demonstrated the activation of YAP caused by ATG7 deletion 74 , which is an important transcription activator in malignant mesothelioma 75 . In our recent research, ATG7's high expression represents a promising prognostic tool for patients with MPM 76 , thus it would be interesting to explore if there is an inverse correlation between ATG7 and YAP in malignant mesothelioma. The leukocyte migration to the site of injury is led by chemokines. The first to be recruited to the site of injury are neutrophils, followed by monocytes, which differentiate into macrophages. Macrophages, once activated, are the main source of growth factors and cytokines which influence the local microenvironment. Mast cells, such as histamine, cytokines, proteases, and lipid mediators also contribute to inflammatory mediators 77 . In previous studies 14, [78][79][80] , it was demonstrated the involvement of cytokines IL-18, IL-1beta, and NF-kappa B in the inflammasome activation process, suggesting that these immune-modulators are involved in the pathogenic mechanisms triggered by asbestos fibers. Acute and chronic inflammations often generate common molecular mediators 70 . The higher inflammation, with greater angiogenesis, causes the fatal outcome of the neoplasm. The release of angiogenic cytokines, including TGF-beta and VEGF, occurs during the angiogenesis process in the malignant mesothelioma progression 81,82 . In particular, VEGF represents the principal angiogenic cytokine involved in this cancer 81,83 , modulating also the development of pleural effusion and ascites through an increase in vascular permeability 84 . VEGF activation plays an essential role in increasing the survival of normal cells exposed to carcinogenic agents. The FE fibers are able to induce functional modifications of parameters with crucial roles in cancer development and progression 9 . The synthesis of VEGF and beta-catenin, two critical steps of epithelial cell activation pathways, are affected by the FE fibers exposure shown by an abnormal cellular status with upregulated cell activities and a risk of neoplastic transformation 85 . Furthermore, the influence of the FE fibers on cell motility has been demonstrated through a dysregulated and altered distribution of actin network 85 . Focal adhesion, which forms mechanical links between cytoskeleton and extracellular matrix (ECM), and adherens junction, result clearly involved in malignant mesothelioma pathogenesis. About that, our recent study, on FE exposure in lung fibroblasts, suggested an ECM remodeling that can give rise to profibrotic cellular phenotypes and tumor microenvironment 86 . Specific signaling pathways have been found to be involved in malignant mesothelioma. Among these, Ras and p53, the commonly mutated genes associated with cancer, are rarely targeted in malignant mesothelioma 87 . Ras has not been found to be mutated in mesothelioma cell lines 88,89 . However, several receptor tyrosine kinase pathways have been shown to be activated in mesothelioma including the epidermal growth factor receptor (EGFR), insulin-like growth factor receptor (IGFR), and c-Met 90-92 , all of which activate Ras signaling. According to these results, several studies have already suggested that the PI3K-Akt pathway is hyperactivated in mesothelioma cell lines 87,93 , resulting in the gain or loss of function of its downstream proteins, 4E-BP1 and pS6, both crucial to the regulation of protein synthesis 94 . But, the prognostic role of the PI3K pathway in MPM is not yet defined 95 .
From papers analyzed by NetME to infer the knowledge network, it is clear that many miRNAs have the ability to positively or negatively influence various pathways involved in apoptosis, necrosis, cell cycle, and cell growth inhibition in malignant mesothelioma, in accordance with our experimental results. Among these,… Furthermore, several pathways, emerging from the literature involved in asbestos exposure, are in common with our data obtained from the expression profiles of miRNAs of healthy mesothelium and MPM in vitro both exposed and not exposed to FE fibers. Among the signaling pathways emerging from our in vitro analyses, we found that some of them are in common with the results processed by NetME. These include e.g. the NF-kappa B, p53, …. already discussed above. Surely, the computational analysis performed to establish the functional role of miRNAs in MPM pathogenesis has shown that these miRNAs can target genes playing a potentially key role in inflammatory processes, angiogenesis, and tumor cell development. Furthermore, miRNAs reach their biological impact by targeting multiple genes with similar biological effects 96 thus, finding potential biomarkers is a very difficult aim that provides a lot of research work. However, surprising results emerged from the expression analysis of tRNA-derived ncRNAs, which have not yet aroused interest in the literature regarding their use as new potential biomarkers in malignant mesothelioma.
After this first preliminary work, the analysis will be followed by the validation of the most significantly deregulated miRNAs by RT-qPCR in an independent sample set. Subsequently, this screening of microRNAs and tRNA-derived ncRNAs will allow us to validate the most interesting and promising results as potential biomarkers for MPM. Our goal is the validation of the most promising results in patients chronically exposed to FE using the liquid biopsy, to provide a minimally invasive screening method for the secondary prevention www.nature.com/scientificreports/ of MPM. Early detection of circulating tumor biomarkers represents one of the most promising strategies to enhance the survival of cancer patients by increasing treatment efficiency 61,97 . Besides this large amount of data, further studies will be designed for the selection of the most significant miRNAs and tRNA-derived ncRNAs to test and validate their diagnostic potential in high-risk individuals. In vitro treatments for RNA-Seq transcriptome profiling. FE fibers were collected from the Biancavilla quarry (Sicily, Italy). These were sterilized under UV light for 10 min, suspended in RPMI 1640 medium, and sonicated through Omni-Ruptor 4000 Ultrasonic Homogenizer (OMNI INternational Inc.; Kennesaw, GA, USA) for 10 min 10 . The stock solution was then diluted appropriately to obtain the different concentrations for in vitro treatments.

Methods
MeT-5A and JU77 were plated onto 100 × 20 mm Petri Dishes (Eppendorf; Hamburg, Germany) at the density of 1 × 10 6 cells and 8.5 × 10 5 cells, respectively. After 24 h of incubation, the medium of both cell lines has been removed and replaced with FE fibers solutions to final concentrations of 50 and 10 ug/ml. MeT-5A and JU77 cells grown in normal medium were used as controls. After 48 h from FE fibers exposure, pellets have been collected in duplicate.
After elimination of the supernatant, cells were harvested on ice by scraping in cold Dulbecco's Phosphate-Buffered Saline (DPBS) (Corning; Manassas, VA, USA). Cells are then centrifuged at 0.2 × g for 5 min at 4 °C and suspended in 1 ml cold DPBS. The cell solution was transferred to Eppendorf tubes. Cells were centrifuged at 0.8 × g for 5 min at 4 °C and supernatant has been removed 10 . The samples have been stored to-80 °C until RNA isolation.

RNA isolation.
Total RNA containing small non-coding RNA was extracted from the cell lines using miRNeasy Mini Kit (QIAGEN; Venlo, Netherlands) according to the manufacturer's recommended protocols (miRNeasy Mini Handbook 11/2020). The RNA was directly quantified by the absorbance ratio at λ = 260/280 nm through NanoDrop (ND 1000) UV-Vis spectrophotometer 10,98 . All samples were diluted at the final concentration of 50 ng/µl for the subsequent analysis.
Small RNA-Seq. QIAseq miRNA library kit (QIAGEN, Hilden, Germany) has been used for small RNA-Seq library preparation following the manufacturer's instructions. RNA samples were quantified and quality tested by Agilent 2100 Bioanalyzer RNA assay (Agilent Technologies, Santa Clara, CA). Libraries were then checked with both Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA) and Agilent Bioanalyzer DNA assay or Caliper (PerkinElmer, Waltham, MA). Finally, libraries were prepared for sequencing and sequenced on singleend 150 bp mode on NovaSeq 6000 (Illumina, San Diego, CA).
Small RNA-Seq analysis. Low-quality reads and adapters were trimmed using Trim Galore, which is a wrapper of FASTQC 99 and Cutadapt 100 . After that, for miRNA analysis, trimmed reads were aligned onto the reference human genome (HG38 version) using HISAT2 101 . The generated SAM files were first converted into BAM files using Samtools 102 , and second, mapped reads were quantified by featureCounts (parameters: -d 14 -primary) using the GTF annotation file retrieved from miRBase 103 . All the above-mentioned steps for the miRNA analysis were performed using RNAdetector 104 software. Concerning the tRNA-derived ncRNA analysis, we first generated an indexed transcriptome using Bowtie 2 105 containing tsRNA and 5' leader pre-tRNA sequences retrieved from tRFexplorer 106 . After that, trimmed reads were mapped on the indexed transcriptome using Bowtie 2 in order to filter them from the reads that may map on tRFs. Filtered mapped reads in SAM format were then converted into BAM files by Samtools 102 , and tsRNAs and 5' leader pre-tRNA sequences were quantified by using bedtools (multicov -q 30). Unfiltered reads in FASTQ format that did not map on tsRNA and 5' leader pre-tRNA sequences were given as input to MINTmap 107 for identifying and quantifying reads associated with tRFs. After that, the raw counts obtained from the miRNA and tRNA-derived ncRNAs analyses were all harmonized together in order to have a single raw count matrix that can be used for the differential expression analysis. The raw count values were then normalized to scale the raw library sizes in Trimmed Mean of M values (TMM) by using edgeR 108 and all miRNAs and tRNA-derived ncRNAs whose geometric mean of TMM values across all samples were less than one were removed from the analysis because they were non-expressed or expressed at a very low level. Finally, the filtered count matrix was used for the differential expression analysis using LIMMA 109 . miRNA and tRNA-derived ncRNAs with a |Log2FC|> 0.58 and an adjusted p-value (Benjamini-Hochberg correction) < 0.05 were considered differentially expressed. Finally, the impact of differentially expressed miRNAs on biological pathways was evaluated by using MITHrIL 50  www.nature.com/scientificreports/ adjusted p-values without a Log2FC cutoff in order to evaluate also the impact of slightly differentially expressed miRNAs on biological pathways.

Determination of dose-response curves and IC50. For the determination of dose-response curves,
MeT-5A were plated onto 96-well plates at the density of 6 × 10 3 cells/50 μl while JU77 were plated at the density of 4 × 10 3 cells/50 μl. After 24 h of incubation, 50 μl of FE fibers solutions were added to the cell cultures in amounts corresponding to final concentrations of 200, 100, 50, 25, 12.5, 6.25, 3.12, 1.56, 0.78 μg/ml. Both cell lines grown in FE-free medium were used as controls. At each time point (6,24,48,72 h of FE exposure) in cell culture, 10% MTT in DPBS has been added to each well. After 4 h of incubation, the lysis solution has been added to each well. The optical density was measured with an absorbance microplate reader at λ = 620 nm. For each sample, 3 replicates were performed. Cell viability was calculated as the percentage of viable cells exposed to FE fibers vs.. control cells no exposed as follows: Cell viability (%) = [OD (Treatment) − OD (Blank)]/[OD (Control) − OD (Blank)] × 100. IC50 values have been calculated through the following equation: [Inhibitor] vs. normalized response-Variable slope.
Results have been analyzed using PRISM GraphPad v. 7.00 and data were represented as the mean ± SD. An unpaired Student's t-test was used to compare data between the two groups. A value of p < 0.05 was considered statistically significant.

Data availability
Raw sequencing data will be deposited on NCBI SRA while processed data will be deposited on GEO". Raw data have been uploaded on SRA, BioProject ID: PRJNA812748.